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Abstract 

In this paper, we consider the short time classical solution to a simplified 
hydrodynamic flow modeling incompressible, nematic liquid crystal materials in 
dimension three. We establish a criterion for possible breakdown of such solu- 
tions at a finite time. More precisely, if (u, d) is smooth up to time T provided 
that J Q T || V x u(t, -)bMO(R3) + ||Vd(t, OHt^s)^ < oo. 
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1 Introduction 

We consider the following hydrodynamic system modeling the flow of liquid crystal 
materials in dimension three (see [2,3,8,10]) and references therein): 

u t + u- V-u - uAu + VP = -Vd ■ Ad, (1.1a) 

d t d + u- Vd = Ad+ \Vd\ 2 d, {1.1b) 

V-u = \d\ = l, (1.1c) 

for (t, x) G [0, +oo) x R 3 . Here u : R 3 — > M. 3 represents the velocity field of the 
incompressible viscous fluid, v > is the kinematic viscosity, P : R 3 — > R represents 
the pressure function, d : R 3 — > S 2 represents the macroscopic average of the nematic 
liquid crystal orientation field, V- and A denote the divergence operator and the laplace 
operator respectively. We are interested in the Cauchy problem (1.1) with the initial 
value 

(u(x, 0), d(x, 0)) = (u (x),d (x)) (1.2) 
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satisfying the following compatibility condition: 



V • uq(x) = 0, |(io(x)| = 1, lim do(x) = a G § 2 , 



(1.3) 



where a is a given unit vector. 

The above system is a simplified version of the Ericksen-Leslie model, which reduces 
to the Ossen-Frank model in the static case, for the hydrodynamics of nematic liquid 
crystals developed during the period of 1958 through 1968 [2,3,8]. It is a macroscopic 
continuum description of the time evolution of the materials under the influence of 
both the flow field u(x,t), the macroscopic description of the microscopic orientation 
configurations d(x,t) of rod-like liquid crystals. Roughly speaking, the system (1.1) is 
a coupling between the non-homogeneous Navier-Stokes equation and the transported 
flow harmonic maps. Due to the physical importance and mathematical challenges, the 
study on nematic liquid crystals has attracted many physicists and mathematicians. 
The mathematical analysis of the liquid crystal flows was initiated by Lin [9], Lin and 
Liu in [10,11]. For any bounded smooth domain in M 2 , Lin , Lin and Wang [12] have 
proved the global existence of Leray-Hopf type weak solutions to system (1.1) which 
are smooth everywhere except on finitely many time slices (see [4] for the whole space). 
The uniqueness of weak solutions in two dimension was studied by [13,15]. Recently, 
Hong and Xin [5] studied the global existence for general Ericksen-Leslie system in 
dimension two. However, the global existence of weak solutions to the incompressible 
nematic liquid crystal flow equation (1.1) in three dimension with large initial data is 
still an outstanding open question. 

The local well-posedness of the Cauchy problem of system (1.1) is rather standard 
(see [4,6,12]). At present, there is no global-in-time existence theory for classical solu- 
tions to system (1.1). In this paper, we will consider the short time classical solution 
to (1.1)-(1.2) and some criterion that characterizes the first finite singular time. Mo- 
tivated by the famous work [1], Huang and Wang [6] have obtained a BKM type 
blow-up criterion. However, the techniques involved in this paper are much different 
from [6], which we believe that the result may have its own interest. Our main results 
are formulated as the following theorem: 

Theorem 1.1. For u G H S (R 3 ) with V • u = and d - a G H S+1 (M. 3 ) with 
|g?o| — 1 for s > 3. Suppose that (u,d) is a smooth solution to the system (1.1)-(1.2), 
then for given T > 0, (u, d) is smooth up to time T provided that 



Notations. We denote by L p , W m ' p the usual Lebesgue and Sobolev spaces on M 3 
and H m = W m ' 2 , with norms || • ||^ P , || • \\w™-p and || • \\ H m respectively. For the sake of 
conciseness, we do not distinguish functional space when scalar-valued or vector- valued 
functions are involved. We denote V = d x = (di, <9 2 , <9 3 ), where di = d Xi) Vj = di and 




(1.4) 
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put d l x f = V' / = V(V Z_1 /). We assume C be a positive generic constant throughout 
this paper that may vary at different places and the integration domain R 3 will be al- 
ways omitted without any ambiguity. Finally, (•, •) denotes the inner-product in L 2 (M 3 ). 

Remark 1.1. It is standard that the condition (1.3) is preserved by the flow. In 
fact, first notice that the divergence free of the velocity field u can be justified by the 
initial assumption that V • u = 0. Indeed, this can be easily and formally observed 
by take V- to the momentum equation. Moreover, applying the maximum principle to 
the equation for \d\ 2 , one also can easily see that \d\ = 1 under the initial assumption 
that | do | = 1- 



2 Proof of Theorem 1.1 

We prove our theorem in this section. Without loss of generality, we assume v — 1. 
The first bright idea to reduce many complicated computations lies in that we just 
need to do the lowest order and highest order energy estimates for the solutions. This 
is motivated by the following observation: 

ii/ii^<cii(/,vv)ii| 2 , v/G# fc . (2.i) 

This inequality (2.1) can be easily proved by combing Young's inequality and Gagliardo- 
Nirenberg's inequality. 

||V7llL P <Cb)||/||2,||vVlll7 a , VfeH k (2.2) 

where - - | = -a + (- - |)(1 - a) with i < k. 

In order to prove Theorem 1.1, we also need the following logarithmic Sobolev's 
which is proved in [7] and is an improved version of that in [1]. 

Lemma 2.1. For p > n, the following logarithmic Sobolev's embedding theorem 
holds for all divergence free vector fields: 

HV/IUoc < C[l + ||/|U2 + || V x f\\ BMO ln(l + ||/|| w2 ,p)]. (2.3) 
Now we are in a position to prove our Theorem 1.1. 



Proof of Theorem 1.1 First of all, for classical solutions to (1.1)-(1.2), one has 
the following basic energy law: 

IK*, -)H! 2 + \\vd(t, .)\\i 2 + jj(||v«( s , oili* + ||Ad( a , •) + \vdMs, .)\\i 2 )ds 

Let's concentrate on the case s = 3. For each multi-index a with |a| < 3, by 
applying d% to (1.1a) and d^ +1 to (1.1b), multiplying them by d^u, d" +1 d respectively 
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and then integrating them over M 3 , we have 



||II^KVd)||| 2 + ||^(V^Arf)||| 2 
= -(d%(u-Vu),d%u) - (d£(Ad- Vd),<9» 

-(d« +1 (u • Vd), d« +1 d) + (^ +1 (|Vd| 2 d), ^ +1 d) (2.5) 

4 

= E ^|a|,i- 



i=l 



where I\ a \^ are the corresponding terms in the above equation which will be estimated 
as follows. Now for \a\ = 1 in (2.5), integrating by parts and using the divergence free 
condition V • u = 0, we have 

IVI = \(dl(u- v«),^«>| = |<^«. v«,^«>| < ciivtiiu-nvtiii^. 

Combining Cauchy's inequality, Sobolev's inequality and the fact |Vd| 2 = — d • Ad 
(since |d| = 1) gives 

|7 li2 1 = | (d x (Ad ■ Vd) , a» I = I (Ad • Vd, 9» | 
<C||Ad||| 3 + |||V 2 9W ||l 2 

< C||Vd||f 4 ||VAd||f 2 + §||V 2 u||| 2 
<G||Vd||i 2 +|(||VAd||| 2 + ||V 2 M ||i 2 ). 

Similarly, 

|/ 13 |= \(d 2 x u- Vd + d x u- Vd x d,d 2 x d}\ 

< C(\\Vu\\ L ~\\V 2 d\\ 2 L2 + ||Vd||i 2 ) + §(||VAd|| 2 2 + \\V 2 u\\ 2 L2 ), 

\h,4 = \(dl(\Vd\ 2 d),d x Ad)\ < C||Vd||^ 4 + §||VAd|| 2 2 . 
Taking the above estimates in (2.5) for \a\ = 1, we arrive at 

|||(V W ,V 2 d)||i 2 + ||(V 2 M ,V 3 d)||i 2 

< C7[||v«iu-(||v«||i a + ||v 2 d||i 2 ) + ||vd||i 2 ]. [Z - b) 

Next we derive an estimate for ||Vd||| 4 . Differentiating (1.1b) with respect to x, we 
have 

d x d t - Ad x d = d x {\Vd\ 2 d-u- Vd). (2.7) 
We multiply (2.7) by \Vd\ 2 d x d and integrate over M 3 to obtain that 

illivd|ii 4 + |||vd||v 2 d|||| 2 + i||4(|vd| 2 )|ii 2 

= (\Vd\ 2 d x d,d x (\Vd\ 2 d-u- Vd)) 

= \\Vd\\ 6 L6 + {\Vd\ 2 d x d,2VdVd x dd- d x u-Vd)) (2.8) 

< C(||Vd|| 3 3 + ||Vu|Uoo||Vd||£ 4 ) 

< C(||V«||^||Vd||l 4 + ||Vd||i 2 ) + |||V 3 d||i 2 . 



We add (2.8) to (2.6) to obtain 

i(\\(Vu, V 2 d)||i 2 + ||Vd|| 4 4 ) + ||(W V 3 d, |Vd||V 2 d|)||| 2 

< c(\\vu\\ L ~ + iivdHiOdKVti, v 2 d)||i 2 + ||vd|| 4 L4 ), 

which gives that 

\\(vu,v 2 d)(t, .)\\i 2 + \\vd(t, -)|| 4 4 + & ||(vy v 3 d, |vd||v 2 rf|)|ii 2 rf s 

< c7c crJ ?» (l|Vtt|l£ -- , - ||Vd|l '* ),to (||(vti, v 2 d)(*o, Olli" + livd(*b, OHIO- 

Noting the condition (1.4) in Theorem 1.1, thus for any fixed small constant e > 0, 
there exists T*(e) < T such that 

/ (||Vx«|| B Afo + ||Vd||i4)ds<e. (2.10) 

Define the temporal energy functional: 

H(t)= sup (||VM S ,-)||| 2 + ||V 4 ^-)||| 2 ) 

T t <s<t 

By (2.1), (2.3), (2.4), (2.9) and (2.10), we arrive at 

ll(v M ,v 2 d)(t,oili 2 + l|vrf(t,oill 4 + /;il(vvvH|vrf||v 2 rf|)|ii 2 ^ 

< ^e/^CollvxuiisMoMi+UKVd)!!^)^ (2.11) 

< C*e CoElQ ( 1+H W = C*(l + H(t)) Co£ , V T, < t < T, 

where Co depends on || (-u , Vdo) III2 and C* depends on ||(Vu, V 2 d)|| 2 2 + ||Vrf||^ 4 . 

Next for \a\ =3. For / 31 , we need to use the following Moser-type inequality (see 
[14, p. 43]): 

\\D s {fg)\W < cdi^iUooiiwil^ + ll/ll^liv^iuo- (2.12) 

Thus we have 

1-^3 i| = \(dldiv{u (g) u), <9 4 w)| 

<C\\V\u®u)\\l 2 + MvMh 

< C\\u\\U\^h + (2.13) 
<C{||VH|| 2 ||V 4 M ||| 2 ||V M ||| 2 ||V 4 M ||| 2 } 2 + ^||V 4 M ||i 2 
<C||V M ||i 4 2 + |||V^||| 2 . 

For J 3i2 , we apply (2.12) and the fact 

1 



Ad ■ Vd = V • (Vd <g> Vd - ^\Vd\ 2 I) 



to obtain that 



|/3,2| = |(d 2 V-(Vd®Vd-l|Vd| 2 /),<9 4 «>| 



+ li- 



< C||V 3 (Vd(g) Vrf- i|V^| 2 J)||| 2 + |||V 4 ||| 
<C||Vd||i.||V 4 d||i a + |||V*«||| 3 (2.14) 

< C{||V 2 d||l 2 ||V 5 d||| 2 ||V 2 d||| 2 ||V 5 d||| 2 } 2 + |||v 4 «||i 2 
<^l|V 2 d||i 4 +|(||V 4 M ||i 2 + ||V 5 d||i 2 ). 
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(2.15) 



Similar as the proof of (2.14), 7 3;3 , / 3i4 can be bounded as follow: 

i4si =m(u-vd),d*d)\ 
<c\m u -vd)\\i 2 + ±\\v 5 d\\i 2 

< c(\\vd\\u\vM\h + IMIUv 4 *) + M^Mlh 

< C{||V 2 rf||l 2 ||V 5 rf||| 2 ||VM||J 2 ||V 4 M||| 2 
+ ||V W ||| 2 ||VS||i 2 ||V 2 rf||| 2 ||V 5 rf||| 2 } 2 + ^|| V 4 W ||| 2 

< c(\\v>d\\% + l|v«||2) + i(\\^u\\h + llv^lliO, 

and 

=\m\vd\ 2 d),d 5 x d)\ 
<c\mwd\'d)\\i 2 +i 2 \\w 5 d\\i 2 

< C(||V 3 |Vrf| 2 ||i 2 + |||Vd| 2 ||i„||V s d||£ a ) + ^||V 5 d||i 2 

< C(||Vd||UV 4 d||i 2 + ||Ad|U|V 3 d||i 2 ) + ^|| V 5 d||| 2 ( 2 - 16 ) 

< c(\\v*d\\ii + {Iiv 2 d||i 2 ||v 5 d||l 2 ||v 2 d||l 2 ||v 5 d||| 2 } 2 ) + M^ 5d Wh 

<C||V 2 d||i 4 2 + |||V 5 d||i 2 . 

Putting (2.13)-(2.16) into (2.5) for \a\ = 3 and by (2.11), we arrive at 

|(l|V 3 n||| 2 + ||V 4 d||| 2 ) + ||V 4 W ||i 2 + ||V 5 d||i 2 
< C{\\V 2 d\\^ + ||V«||i 4 2 ) < C*(l + H(t)f^ 

for all < t < T. Integrating the above inequality with respect to time from to 
t G [T*, T), we obtain 



(2.17) 



i + l|vM^-)lli 2 + llvMt,-)lli 2 

< 1 + || VMT*, .)||| 2 + || V 4 d(T„ -)||| 2 + C.(l + H(t)f^. 
If we choose e small such that C £ < |, then (2.17) implies 

i + < c(i + 1| v 3 «(r„ -)lli 2 + II v 4 d(T„ oilio, 

for all T* < t < T. Since the right hand side of above inequality is independent of t for 
T* < t < T, we conclude that the above inequality is valid for t = T which means that 
u(T, ■) G iJ 3 (lR 3 ) and d(T, •) — a G iJ 4 (lR 3 ). Thus the proof of Theorem 1.1 is completed. 
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